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I.  INTRODUCTION 


Small  laboratory  devices  such  as  vented  chambers^or  blowout  guns 
have  a  long  history  in  understanding  gun  barrel  wear.  Compared  with 
large-caliber  gun  tests,  such  devices  need  less  propellant,  need  less 
room  for  safe  firing,  and  have  inserts  which  can  be  removed  for  post¬ 
firing  analysis.  An  unanswered  question  is  how  to  relate  wear  in  the 
devices  with  wear  in  large-caliber  guns.  A  recent  .JANNAF  workshop  con¬ 
cluded  that  laboratory  devices  cannot  evaluate  wear-reducing  additives, 2 
and  that  test  devices  in  three  laboratories  have  given  contradictory 
results  with  the  same  propellants. 3 >4 

This  report  attempts  to  relate  erosivity  from  propellants  in  the 
BRL  37mm  blowout  gun  to  wear  expected  in  large  caliber  guns.  Empirical 
formulas5,6  presume  the  wear  is  exponentially  proportional  to  the  maxi¬ 
mum  bore  surface  temperature.  Bore  surface  temperatures  in  the  empira- 
cal  models  are  estimated  with  interior  ballistic  parameters  such  as 
charge  weight,  flame  temperature,  gun  diameter,  and  peak  pressure. ' 


1 

''Hypervelocity  Guns  and  the  Control  of  Gun  Erosion ,  "  Summary  Technical  ■ 
Report  of  Division  1 ,  NDRC,  Volume  1,  Washington3  DC  1946. 

2 

J.A.  Lannon  and  J.R.  Ward3  "Workshop  Report  on  Wear-Reducing  Additives 
and  Their  Performance  in  Guns 3  "Proceedings  of  the  17th  JANNAF  Combustion 
Meeting3  CPIA  Publication  329,  November  1980. 

3 

A.J.  Bracuti3  L.  Bottei3  J.A.  Lannon 3  and  L.H.  Caveny3  "Evaluation  of 
Propellant  Erosivity  with  Vented  Chamber  Apparatus 3  "  1980  JANNAF  Pro¬ 
pulsion  Meeting 3  CPIA  Publication 3  315,  March  1980. 

^ J.R.  Ward3  R.W.  Geene3  A.  Niiler3  A.  Rye3  and  B.B.  Gr oilman,  "Blowout 
Gun  Erosivity  Experiments  with  Double-Base 3  Triple-Base 3  and  Nitramine 
Propellants 3  1980  JANNAF  Combustion  Meeting  CPIA  Publication 3  315 3 
March  1980. 

5j. M.  Frankie  and  L.R.  Kruse 3  "A  Method  for  Estimating  the  Service  Life 

of  a  Gun  or  Howitzer 3"  BRL  Memorandum  Report  No.  1852 ,  June  1967. 

6  (AD  #818348)  . 

C.S.  Smith  and  J.S.  O'Brasky 3  "A  Procedure  for  Gun  Barrel  Life  Estima¬ 
tion,  "  Proceedings  of  the  Tri-Service  Gun  Tube  Wear  and  Erosion  Sympo¬ 
sium,  ADPA,  Dover,  NJ,  March  1977. 

7 

J.  Comer,  Theory  of  the  Interior  Ballistics^  of  Gunsj  John  Wiley  &  Sons , 
Inc.,  NY,  1950. 

Q 

L.W.  Nordheim,  H.  Soodak,  and  G.  Nordheim,  "Thermal  Effects  of  Propellant 
Gases  in  Erosion  Vents  and  Guns,"  NDRC  Armor  and  Ordnance  Report  No. 
A-262,  March  1944. 
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This  report  uses  Nordheim’s  scheme  to  compute  the  peak  surface 
temperatures  during  wear  tests  of  five  propellants . ®  If  the  relative 
wear  of  the  propellants  in  the  blowout  gun  is  the  same  as  in  large  guns, 
then  a  plot  of  natural  logarithm  of  wear  v£  peak  temperature  should  be 
linear,  with  the  slope  identical  to  that  in  guns.  Wear  tests  ip  Reference 
9  incorporated  Niiler's  radioactive  technique^  to  correlate  mass  loss  and 
diameter  change. 

II.  INPUT  PARAMETERS  AND  COMPUTER  CODES 

Three  computer  codes  were  used  to  perform  the  calculations  dis¬ 
cussed  in  the  report.  The  first  program,  IB  3/70,  computes  interior 
ballistics  data  for  guns.  The  second  program,  NOZZLE,  computes  interior 
ballistic  data  for  a  blowout  gun  with  an  erosion  nozzle.  The  third  pro¬ 
gram,  NEWNSN,  computes  heat  transfer  using  files  of  IB  3/70  and  NOZZLE. 
Listings  for  each  program  are  contained  in  the  Appendix. 

The  assumed  heat  input,  as  described  fully  in  Reference  8,  was 
taken  to  be 


Q  =  JgXCppUAT  ,  (1) 

2 

where  q  =  heat  flux,  J/m  *s, 

X  =  friction  factor, 

C  =  specific  heat,  J/kg*K, 

3 

p  =  propellant  gas  density,  kg/m  3 
U  =  gas  velocity,  m/s,  and 

AT  =  temperature  difference  between  propellant  gas  and  wall,  K. 

Following  Nordheim's  method,  a  friction  factor  of  1/253  was  used 
for  the  17.3  mm  diameter  nozzle. 

Computations  were  done  for  the  five  propellants  evaluated  in  Re¬ 
ference  9.  The  compositions  are  listed  in  Tables  1  and  2. 


9  R.W.  Geene ,  J.R.  Ward 3  T.L.  Brosseau ,  A.  Wilier ,  R.  Birhnire,  J.J. 
Rocchio,  "Erosivity  of  a  Nitramine  Propellant , "  BRL  Technical  Report 
No.  020943  August  1978.  (AD  #A060590) 

10S.E.  Caldwell  and  A.  Niileri  " The  Measurement  of  Wear  from  Steel  Using 
the  Radioactive  bbCoy"  BRL  Report  No.  1923 3  September  1976.  (AD  #A030262) 
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TABLE  1.  COMPOSITIONS  AND  GRAIN  DIMENSIONS  OF  M30  AND  HFP  PROPELLANTS 


M30  HFP 

PPL-A-6372  PPL-A-6380 


Nitrocellulose  (12.6%N) 

28.0% 

29.3% 

Nitroglycerin 

22.5 

22.7 

Nitroguanidine 

47.7 

5.0 

RDX 

36.5 

Dioctylphthalate 

5.0 

Ethyl  Centralite 

1.5 

1.5 

Cryolite 

0.3 

Total  Volatiles  (Residual) 

0.2 

0.3 

Grain  Length,  mm 

7.78 

10.58 

Grain  Diameter,  mm 

1.59 

2.37 

Grain  Perf.  Diameter,  mm 

0.46 

0.77 

Grain  Web,  mm 

0.56 

0.80 

Grain  Geometry 

SP 

SP 

TABLE  2.  COMPOSITIONS  AND 

GRAIN  DIMENSIONS 

OF  M5 ,  M8 ,  AND  Ml 

PROPELLANTS 

M5_ 

M8 

Ml 

Nitrocellulose  (13.25%N) 

81.95% 

52.15% 

85.00% 

Nitroglycerin 

15.00 

43.00 

Ethyl  Centralite 

0.60 

0.60 

Barium  Nitrate 

1.40 

Potassium  Nitrate 

0.75 

1.25 

Diethylphtalate 

3.00 

Dinitrotoluene 

10.00 

Dibutylphthalate 

5.00 

Diphenylamine,  Added 

1.00 

Ethyl  Alcohol,  Residual 

2.30 

0.40 

0.75 

Water,  Residual 

0.70 

0.50 

Graphite 

0.30 

Grain  Length,  mm 

10.58 

25.4 

8.26 

Grain  Diameter,  mm 

3.92 

12.7 

3.68 

Grain  Perf.  Diameter,  mm 

0.41 

0.37 

Grain  Web,  mm 

0.69 

0.56 

0.64 

Grain  Geometry 

7  Perf 

Strip 

7  Perf 
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Propellant  parameters  required  for  the  computations  are  listed  in 
Table  3. 

The  term,  B/W,  is  used  to  relate  burning  rate,  r,  with  chamber  pressure,  P, 
and  propellant  web,  W,  by 

r  =  (P  +  44.8) (J)  •  W  .  (2) 

Propellant  webs  are  listed  in  Tables  1  and  2. 

3 

For  the  blowout  gun,  chamber  volume  and  nozzle  area  were  315  cm 
and  2.36  cm2,  respectively. 


III.  RESULTS  AND  DISCUSSION 

The  experimental  results  from  Reference  9  are  listed  in  Table  4. 

Table  5  gives  calculated  time  to  maximum  pressure,  peak  nozzle  sur¬ 
face  temperature,  heat  input,  and  time  from  nozzle  rupture  to  peak  nozzle 
surface  temperature. 

A  linear  least  squares  fit  tested  the  exponential  dependence  of  wear 
with  peak  surface  temperature. 

£n  (w)  =  a  +  bT  ,  (3) 


where  w  =  wear,  y/shotj 

T  =  peak  bore  surface  temperature,  K,  and 
a,b  =  constants. 

The  best  fit  values  of  a  and  b  were  -9.7  ±  2  and  0.0094  ±  0.002, 
respectively,  with  error  representing  one  standard  deviation.  Table  6 
compares  experimental  wear  and  wear  computed  with  Eq.  (3).  Eq.  (3) 
clearly  underestimates  wear  for  Ml  with  a  peak  temperature  of  773  K. 
Frankle-Kruse  also  noted  their  expression  was  invalid  below  900  K; 
Smith-0' Brasky  report  a  similar  threshold  of  750  K. 

A  check  on  how  wear  varies  with  peak  bore  surface  temperatures  in 
guns  can  be  made  with  bore  surface  temperature  computations  by  Vassallo 
and  coworkers  at  Cal span  Corporation! 1  in  the  M185  cannon.  The  Cal span 
workers  compute  peak  bore  surface  temperatures  from  measured  total  heat 
input. 12  Recent  work  in  a  shock-tube  gun  with  different  metals  showed 
temperature  calculations  starting  with  measured  total  heat  input  could 
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TABLE  3.  PROPELLANT  PARAMETERS* 
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TABLE  4.  EXPERIMENTAL  DATA  FROM  BLOWOUT  GUN  EXPERIMENTS 
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TABLE  5.  COMPUTED  HEAT  TRANSFER  RESULTS 
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TABLE  6 


COMPARISON  BETWEEN  EXPERIMENTAL  WEAR  AND  WEAR  COMPUTED 
WITH  EQUATION  (3) 


Propellant  Peak  Temp. ,  K  Wear  exp't,  y/Shot  Wear  Cal’d,  n/Shot 


Ml 

Ml 


773 

875 


0.2 

0.1 


0.0005 

0.2 


M30 

966 

0.5 

M30 

1,130 

0.6 

M30 

1,228 

3.8 

0.5 

2.5 

6.2 


HFP 

968 

0.5 

HFP 

1,150 

1.1 

HFP 

1,235 

6.8 

0.5 

3.0 

6.6 


M5 

982 

0,8 

M5 

1,155 

4.1 

M5 

1,250 

18.5 

0.6 

3.1 

8.3 


M8 

1,107 

M8 

1,210 

MS 

1,407 

2.8 

9.6 

48.6 


2.0 

5.0 

32.1 
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predict  correctly  whether  the  melting  point  of  the  metal  was  exceeded. 
Starting  conditions  in  the  shock-tube  gun  were  varied  to  produce  tempera¬ 
ture  bracketing  the  metals'  melting  temperatures. 

Table  7  shows  bore  surface  temperatures  computed  for  three  155mm 
propelling  charges  which  have  no  wear-reducing  additive.  Wear  data  is 
available  for  the  M119  and  the  XM201E2  charges, 14 >15  Wear  sensor  data  in 
the  M185  cannon  showed  the  XM208  minus  its  additiye  was  five  times  as  ero¬ 
sive  as  the  XM201E2.  Wear  of  the  XM208  was  presumed  to  be  five  times 
greater  than  the  measured  wear  of  the  XM201E2  charge.  A  least-squares 
fit  of  the  three  points  in  Table  8  to  Eq.  (3)  produced  values  of  a  and  b 
of  -7.4  and  0.0076K-!,  respectively. 

TABLE  7.  BORE  SURFACE  TEMPERATURES  AND  WEAR  FOR  155mm  CHARGES  WITHOUT 
ADDITIVE 


Charge  Peak  Temp,  K  Wear,  p/shot  Ref. 


M119 

961 

0,9 

15 

XM201E2 

1,100 

2.6 

16 

XM208  (no  additive) 

1,306 

12.5 

12 

Figure  1  plots  wear  vs  peak  bore  temperature  for  the  nozzle  and  the 
155mm  gun.  Figure  1  suggests  relative  wear  among  propellants  will  be 
overestimated  compared  with  wear  in  guns  by  yen  percent.  Nonetheless, 
the  agreement  seems  close  enough  to  suggest  relative  propellant  erosivity 
can  be  inferred  in  nozzles. 


11 

F.A.  Vassallo,  "An  Evaluation  of  Heat  Transfer*  and  Erosion  in  the  155 
mm  M185  Cannon 3"  Calspan  Technical  Report  Ho.  VL-5Z37-B-1,  July  1976. 

12 

F.A.  Vassallo,  "Mathematical  Models  and  Computer  Routines  Used  in  Ev¬ 
aluation  of  Caseless  Ammunition  Heat  Transfer , "  Calspan  Technical  Re¬ 
port  No.  GM-2948-2-1,  June  1971. 

1 3 

F.A.  Vassallo  and  W.R.  Brown,  "Shock  Tube  Gun  Melting  Erosion  Study 3" 
BRL  Contractor  Report  Ho.  406 3  August  1979.  (AD  #A076219) 

14 

J.J.  Reed  and  J.P.  Cherry 3  "Service  Tests  of  155mm  Howitzer 3  Self- 
Propelled  Equipped  with  XM185  Tube3  "  Field  Artillery  Board  Report. 

Ft.  Sill 3  0K3  January  1970. 

1  5 

J.A,  Bemaree,  "155mm  Ml 8 5  Tube  Wear  Tests  of  Charge 3  Propelling  XM2013 " 
Interim  Report 3  JPG-76-601 3  June  1976. 
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WEAR,  fx  /SHOT 


50.0 


10.0 


5.0 


1.0 
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0.1  I - 1 - 1 - 1 - 1 - 1 — 

0.8  0.9  1.0  1.1  1.2  1.3 

PEAK  BORE  SURFACE  TEMP,  K,  *103 

Figure  1.  Wear  vs_  Peak  Bore  Surface  Temperature 
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Some  other  observations  from  the  calculations  of  bore  surface  tem¬ 
peratures  can  be  made  from  comparisons  of  bore  surface  temperature  vs. 
time  for  the  37mm  blowout  gun  and  a  3"/70  Naval  gun  depicted  in  Figures 
2-4.  Table  8  summarizes  pertinent  heat  transfer  data  for  the  three 
figures . 


TABLE  8.  CALCULATED  HEAT  INPUT  AND  SURFACE  TEMPERATURES  FOR  3"/70 
AND  37mm  BLOWOUT  GUNS 


Gun 

Nozzle 
Diameter , mm 

1/A 

Propellant 

Charge 

Mass,g 

Flame 

Temp,K 

Surface 
Temp , K 

Heat  In-^ 
put, J/mm 

37mm 

17.3 

253 

HFP 

89.8 

3,255 

1,235 

0.78 

3"/70 

- 

280 

Picrite 

3991. 

2,065 

1,224 

1.13 

37mm 

13.3 

239 

Ml 

159 

2,417 

1,254 

1.21 

The  first  point  of  interest  is  the  short  heating  time  of  the  blow¬ 
out  gun  in  Figure  2  vs.  the  3"/70  Naval  gun.  As  a  result  of  the  longer 
heating  time,  the  3"/70  gun  reaches  the  same  peak  bore  surface  tempera¬ 
ture  with  a  much  cooler  propellant.  The  effect  of  the  longer  heating 
time  is  also  reflected  in  the  larger  total  heat  input  for  the  gun.  The 
figures  also  point  out  why  propellants  with  flame  temperatures  near 
3800K  would  be  needed  to  get  the  wear/shot  characteristic  of  tank  guns 
with  the  17.3  mm  nozzle.  Figure  4  shows  that  Ml  propellant  produces  an 
equivalent  peak  temperature  as  HFP  if  charge  mass  is  increased  and  nozzle 
diameter  is  decreased.  A  new  chamber  would  be  required  in  the  37mm  gun 
to  be  able  to  load  160g  of  propellant. 


IV.  CONCLUSIONS 

1.  Bore  surface  temperature  vs.  time  were  computed  for  17.3  mm  diameter 
nozzles  in  the  BRL  37mm  blowout  gun  in  order  to  test  dependence  of  wear 
vs.  peak  bore  surface  temperature. 

2.  The  natural  logarithm  of  wear  from  the  nozzle  was  shown  to  be  linearly 
dependent  on  peak  bore  surface  temperature  which  agrees  with  empirical 
wear  models  and  with  calculations  of  peak  bore  surface  temperatures  by 
Calspan  Corp.  for  the  155mm  M185  cannon.  It  was  shown  the  slope  is  10% 
steeper  in  the  nozzle  suggesting  relative  propellant  erosivity  with  the 
17.3mm  diameter  nozzle  may  be  exaggerated  slightly  compared  to  large 
caliber  guns. 
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Figure  2.  Temperature  vs.  Time  In  37mm  Blowout  Gun  (17.3mm 
Diameter  NozzTe-HFP  Propellant). 
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TEMPERATURE  (DEG.  K ) 

400  G00  900  1000  12O0  14D0 


Figure  3.  Temperature  vs.  Time  In  3"/70  Gun  With  Plcrite  Propellant. 
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TEMPERATURE  (DEG.  K) 

400  S00  600  100G  1200  1400 


Figure  4.  Temperature  vs.  Time  in  37mm  Blowout  Gun  (13.3mm 
Diameter  NozITe-Ml  Propellant). 
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APPENPIX 


COMPUTER  CODRS 


COMPUTER  CODES 


The  first  program,  "IB3/70",  computes  interior  ballistics  data 
for  guns . 

The  second  program,  "NOZZLE",  computes  interior  ballistics  data 
for  a  blowout  gun  with  an  erosion  nozzle. 

The  third  program,  "newnsn"  computes  heat  transfer  and  conduction 
data.  The  program  "newnsn"  uses  information  stored  on  files  !'DATA" 
and  "DATA1"  by  "IB3/70"  or  "NOZZLE". 
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"IB3/70" 


Input  data  begins  in  line  2860.  Data  to  be  stored  in  each  line 


is  as  follows: 

Line 

Data 

2860 

Projectile  weight  (lbs) 

2870 

Projectile  travel  (inches),  chamber  volume  (cu.in.) 

2880 

Pidduek-Kent  constant,  integration  step  size  (sec) 

2890 

Estimated  projectile  velocity  (ft/sec),  bore  diameter (in) , 
bore  area  (sq.in.) 

2900 

Shot  start  pressure  (psi) 

2910 

Igniter  weight  (pounds) 

2920 

Igniter  impetus  (in  lbf/lbm),  ratio  of  specific  heats, 
flame  temperature  (°K),  covolume  (cu.in. /lb) 

2930 

Propellant  weight  (lbs) 

2940 

Propellant  force  (  in  lbf/lbm) ,  ratio  of  specific  heats, 
flame  temperature  (°K),  covolume  (cu.in. /lb),  density (lb/cu. in 

2950 

Propellant  grain  diameter  (in),  length  (in),  perforation  dia¬ 
meter  (in) ,  and  number  of  perfs 

2960 

Propellant  identification  (an  alphanumeric  string) 

2970 

PI,  B (where  the  burning  rate  (in/sec)  is  B* (space  mean 
pressure  in  psi  +  PI)),  number  of  integration  steps  per 
printout,  expected  peak  chamber  pressure  (psi),  expected 
peak  acceleration  (g's). 

2980 

Number  of  "resistive  pressures'1  listed,  followed  by  the 
appropriate  number  of  pairs  of  travel  distance (in)  and 
resistive  pressure  (lb/sq.in.) 

2990-3020 

These  lines  are  used  by  the  program  to  aid  in  calculating 
the  fraction  of  the  7  perf  grain  burnt  after  splintering. 

3025 

Use  this  line  to  provide  data  if  called  for  by  the  input 
statements 

Data  requested  by  input  statements: 

"Number  of  parameters  to  be  varied":  Enter  a  number  (0  to  5). 

"Number  of  proj  wt,  cham  vol,  pmax,  travel  distance,  charge  wt":  The 
number  of  non-zero  enteries  should  be  as  indicated  in  the  immediately 
previous  input.  The  values  for  the  parameters  indicated  should  have  been 
listed  in  line  3025  of  the  program. 

"Option  No  2" 

Enter  a  number  (0  to  7).  The  program  will  then  perform  as  indicated 
below: 
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NUMBER  ENTERED 


ACTION 


0  Compute  IB  trajector  for  input  data  (i.e.,  data 

in  lines  2860-2980) 

1  Change  B  to  match  indicated  peak  pressure 

2  Change  charge  weight  to  match  indicated  pressure 

3  Change  B  to  match  indicated  acceleration 

4  Change  charge  weight  to  match  indicated  acceleration 

5  Change  charge  weight  to  match  indicated  velocity 

6  Print  out  a  blurb  of  the  above 

7  Same  as  0,  but  program  writes  data  on  files  DATA 
DATA1  for  future  use. 

Note  that  if  the  first  input  requested  was  not  zero  multiple  runs  will  be 
performed.  This  is  not  useful  if  "Option  number  2"  is  7. 
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10. 

OPTION  ERSE  1 

20  O' 

■'ERL  HP 

30 

DIM  C < 50 > 

40 

DIM  P* 1 10] 

50 

DIM  T  <  2 , 1 0  > , F  <  2 , 1 

2  )  ,  R  *s  6  „  3  > 

68 

Ehe-at^O 

7Q 

LET  26=0 

SO 

D E F  F N fi <  Y 9 >  =  IN T  < V 9*1 0 6 0 > s 1 0 0 0 

90 

F:  E  A  D  G 1 ,  G  2 ,  G  3 ,  G  5 , 

G  6  •,  G  7 ,  i  j  3 ,  G  4  , 

100 

LET  E7=G?*12 

110 

READ  I  1 , 12,  13, 14, 

15 

120 

READ  Pi , P 2 , P 3 , P 4 , 

P5,P6,P7,P8, 

130 

READ  Ft 

140 

READ  H , B , M 2 , M  3 , A  9 

150 

IF  P0O7  THEN  190 

160 

LET  H  1  =  < P 7 - 3 * P 9 > 

4 

176 

LET  W2=  <  < P7-W1 >  A2 

■"  7  *  <  p  9  +  W 1  )  2 

ISO 

GOTO  200 

190 

LET  W1 =1090 

2  0  0 

READ  K 

216  FOR  1=1  TO  K 

220  READ  TCI,  I  >,  T.<2,  I  > 

236  NEXT  I 
246  FOR  1=1  TO  12 
250  READ  F(1,I),FC2,D 
260  NEXT  I 

270  PRINT  " r. UMBER  OF  PARAMETERS  TO  PE  VARIED" 

280  INPUT  D 

290  IF  D=0  THEN  516  * 

3O0  PRINT  "NO  OF  FRO-J  NT  ,  CHAM  VOL ,  PM  AX,  TRAVEL  III  ST,  CHARGE  NT" 

310  I NPIJT  D 1  ,  D2 ,  D3 ,  D4 ,  D5 

320  IF  Dl=6  THEN  366 

330  FOR  N 1 = 1  TO  D1 

346  READ  L<N1>. 

350  NEXT  N1 
360  IF  D 2  =  O  THEN  408 
378  F 0 R  N 2=1  TO  D  2 
336  READ  MCN2) 

390  NEXT  N2  ; 

406  IF  D3=0  THEN  470 
410  FOR  N3= 1  TO  D3 
420  READ  Q < N 3  > 

438  NEXT  N3 

440  FOR  N 4  = 1  TO  D4 

450  READ  W < N4 ) 

460  NEXT  N4 
470  IF  D 5  =  0  THEN  510 
480  FOR  N 5  =  Q  TO  D5 
496  READ  Y<N5> 

500  NEXT  N5 

510  PRINT  "OPTION  NO  2  (0OX>»5>  ;  INPUT  6  fOR  EXLANAT  ION" 

520  INPUT  Cl 

530  IF  C 1 < 6  THEN  760 

540  IF  Cl =6  THEN  620 

559  C1=0 

560  CREATE  "DATA" ,126 
570  CREATE  "DATA1" ,3 
586  ASSIGN  #1  TO  "DATA" 

581  BUFFER  #1 

598  ASSIGN  #2  TO  "DATAl" 

6O0  Eh eat = 1 
616  GOTO  930 
620  PRINT 
630  PRINT 
640  PRINT 

650  PRINT  "TO  MATCH" ; TAB <20 ) ; "BY  ADJUSTING"  ;  TAB<41  :  11  INPUT  " 
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660 

670 

680 

690 

7O0 

710 

720 

730 

740 

750 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

INPUT 


“NOTHING"; TAB (43); “Q" 

"PMRX" ; TAB < 20) ; “BETA" ; TflB<43> ; “1" 

" P M R X “ ; TAB<20> ; “CHARGE  NT" ; TAB<43> ; "2" 

“ ACCEL “ ; TAB  < 20 > ; “BETA" ; TAB  < 43  > ; " 3“ 
"ACCEL" ; TAB<20>; "CHARGE  NT " ; TAB < 43 ) ; "4" 


,'ELCC  I ' 


TAB < 20) ; “CHARGE' NT" ; T  A  B  <  4  3 > j 


Cl 


760 

IF  D 

=0  THEN  930 

770 

IF  D1=0  THEN  300 

780 

FOR 

N 1 = 1  TO  HI 

790 

LET, 

G1=L<N1 > 

800 

IF  ri 

2=0  THEN  830 

8  1  O 

FOR 

H2= 1  TO  D2 

820 

LET 

G  3  =  M  <  N  2  > 

830 

IF  Ii 

3 = 0  THE  N  8 6  0 

840 

FOR 

N3= 1  TO  D3 

850 

LET 

M  3  =  Q  (  N  3  > 

860 

IF  D  4  =  0  THEN  890 

870 

FOR 

N4= 1  TO  B4 

8  y  y 

LET 

G2  =  N  <  N4  > 

890 

IF  Ii 

5=0  THEN  930 

900 

FOR 

N5= 1  TO  D5 

910 

LET 

P 1 =Y<N5 > 

9  2  0 

NEXT 

'  N5 

930 

LET 

J 1 = I 1*12/(13-1 > 

940 

LET 

J  2  =  J  1  /  I  4 

950 

LET 

J3= I 1*12/14 

9  6  O 

LET 

J5=I1*I5 

970 

LET 

Q 1  =  P  1  *  P  2  /  <  P  3  - 1  ) 

980 

LET 

Q  2  =  Q 1  P  4 

9  9  0 

LET 

Q3=P1*P2/R4 

1  000 

LET 

Q  4  =  P 1  /  F*  6 

1010 

LET 

Q  5  =  P 1  *  P  5 

1020 

LET 

Q6=I 1+P1 

1030 

Q7  =  F 

>  8  *  <  P  7  A  2  -  P  0  *  P  9  •-  2  > 

1040 

LET 

Q  8 = G  6  /  2 

1  050 

LET 

0 1  =  •  3 y  * i j 8 1 . 5*  <  G2+G 3  /  G 4  )  * 

i  0  6  0 

LET 

0 1 =0 1  /  <  1  + . 6*G8*2 . 175/P  1 * . 

1070 

LET 

H 1 =0 1  +  ( G 1 +Q6/G5 > / 7  7  2 

1  080 

LET 

0 2=1 / <  P  3 -  1  > 

1090 

LET 

H  2  =  1  +  Q  6  /  (  G  1  *  G  5  > 

1100 

LET 

03=  1  / (.  < 2*02+3 >/G5+2* < 02+ 1 

1110 

LET 

H3=a-03>M-02~1  > 

1  120 

LET 

H  4  =  G  4  *  3  8  6  /  G  1 

1130 

IF  M2 >=150  THEN  1160 

<  p  4  -  3  0  0  > 

8  3  7  5  >  /  (.  1 4  4  *  G  7  *  G  7  > 


1140  PRINT  “TIME  BIST 

1150  PRINT  "FR.  ACCEL 

1160  LET  T=1000 
1170  LET  T  2  =  T  <  2 ,  1  ) 

1180  LET  N=0 
1190  LET  Ml =1 

12O0  LET  R ( 1 ,  1 ) =  R ( 1 , 2 ) =  R < 1 , 
1210  LET  X 1  =  R  <  1 , 2 > 

1220  LET  V 1 =R C 1 , 3> 

1230  FOR  1=1  TO  3 
1240  LET  R<2, I >=R< 1 , I > 

1250  NEXT  I 

1260  LET  J  =  2 

127G  GOSUB  2280 

128G  LET  S-S+Q8 

1290  LET  J  =  3 

1300  FOR  1=1  TO  3 

1310  LET  R<2, I)=R<2, I>*G6 


PRES 

BASE" 


VEL 


TEMP 


>=B0=2~T2=S=B3=B4=B5=Z5=26^0 


i 
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1320  LET  R  3  ,  I  )  -R  <  1  ,  1  )  +R  <  2 ,  I  >  /2 

1336  NEXT  I 

1  3  4  0  G  0  S  U  B  2  2  8  8 

1350  LET  J =4 

1360  FOR  1=1  TO  3 

1370  LET  R < 3 ,  I  >  =  R < 3 , I > *G 6 

1380  LET  R<4, I >=R< 1 „ I>+R<3, I)/2 

1390  NEXT  I 

1  4  0  0  G  0  S  U  E  2 2  3  0 

1410  LET  J=5 

1420  LET  S=S+Q8 

1430  FOR  1=1  TO  3 

1440  LET  R<4, I>=R<4, I>*G6 

1450  LET  R <  5  .  I  >  =  R  <  1  ,  I  >  +  R  <  4  ,  I > 

1460  NEXT  I 

14  7 0  GOSU B  2280 

1480  FOR  1=1  TO  3 

1490  LET  R<5, I >=R<5, I >*G6 

1500  LET  R(l, I >  =R < 1 , I )  +  <  R  <  2 , I)  +  R<3, I)*2+2*R<4? J>+R<5f I  >  }/6 
1510  LET  R<6, I >=R<  1 , I > 

1520  NEXT  I 
1530  LET  J=6 

1548  IF  Eh at  THEN  PRINT  #1 ;  V,  X,  Pl*2/<G3+G4*>NPl*U-2>/P6>  ,  T 
1550  IF  M2* I  NT  <  M 1 /M2  >  <  >M 1  THEN  1590 

1 560  PR  1  NT  S* 1 0O0 ; TfiB < 8 > ; FNR < R <1 , 2 ) > ; TRE < 1 8 > ; I NT < 30 > ; TRB  <  27  >  5 
1  57G  PRINT  I  NT  <  R  <1  ?  3  )  /•  1 2  >  ;  TR3  <  34  )  ;  I  NT  <  T  >  ;  TRB  C  45  >  ;  FNR  <  Z  >  5 
1580  PRINT  TRB < 54); I  NT <R 1/386. 09> ; TRB<63> ; ZNT<B2> 

1590  IF  R < 1 , 2 ) > = G 2 - 2 5  THEN  1610 

160O  GOTO  1630 

1610  H  =  N  + 1 

1626  LET  C < N >  =  B2 

1630  IF  R < 1 , 2 > > = G 2  THEN  1730 

1640  LET  M1=M1+1 

165G  LET  V  =  R < 1 , 3 ) 

1660  IF  B 8 <  = E 3  THEN  1210 
1679  LET  R 5  =  R 1 / 3 8 6 . 0 9 
1630  LET  B3=B0 

1690  LET  B 8  =  B 2  ; 

1700  LET  B 4  =  8 
1718  LET  E5  =  RU,2> 

1720  GOTO  1210 

1730  LET  X=<G2-X1 )  /  <  R  < 1 , 2  >~X  1 ) 

1740  LET  V  =  V 1+X*<R(1 , 3 ) - V 1 > 

1750  LET  S 1 =  S - Q 8 

1760  LET  29  =  C < N ~ 1 >  +X*K C < N ) - C < N - 1  >  } 

1770  IF  C 1=0  THEN  I960 
1  7  8 0  B r  an c  h = C 1 

1798  IF  <Branch>=l >  RNB  <Branch<6)  THEN  ON  Branch  GOTO  1808, 1830, I860, 1890, 1920 
1800  IF  RBS < B3-M3 > <200  THEN  I960 
1810  LET  B  =  B/ < 1 - . 85*  <  M3-B3  > / < M3+B3 >  > 

1320  GOTO  1940 

1830  IF  RBS < B3-M3 ) < 200  THEN  I960 
1840  LET  P 1 =P 1 / ( 1 - . 9*  < M3-B3 > / < M3+B3 > > 

1850  GOTO  1940 

I860  IF  RBS<R9-R5><100  THEN  1960 
1870  LET  B  =  E / < 1-1 . 1  * < R9-R5 > s < R9+ R5 > ) 

1880  GOTO  1940 

1890  IF  RBS<R9-R5)< 100  THEN  I960 
19G0  LET  P  1  =p  1  /  1  -  .  9*  <  R9-R5  )  /  < R9  +  R5 )  ) 

1910  GOTO  1940 

1920  IF  RES < G7- V/Q2 ) < 1 0  THEN  I960 

1938  LET  P1=P1/<1-1 .5*(G7-V/12>/(G7+V/12>> 

1940  LET  E7=V 
1950  GOTO  930 
I960  PRINT 

1970  IF  Eheat,=0  THEN  2000  29 


1980 
2000 
2010 
2O20 
203O 
2040 
2O50 
2060 
207O 
2080 
2090 
2100 
2110 
212G 
2130 
2 1  40 
2150 
2160 
2170 
2180 
2  1  90 
2260 
2210 
2220 
2230 
2240 
2256 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
234  0 
2350 
2360 
2370 
2380 
2390 
2490 
2410 
2426 
2430 
2448 
2450 
2460 
2470 
2480 
2490 
25O0 
2510 
2526 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2648 


CALL  F-: j  <#1  ,  # 2  ,  P  1  **''  < 6 3  +  G 2 *64 >  ,  <V>  ,  S, 66, T  ,  P5,  62+63/64 , 68 ,  G3.'G4  > 

IF  M2 >=50O  THEN  2190 

PRINT  "PRO J . NT .  — " ; G 1 ; "LBS . " 

PR  I  NT  “  I  HI  T  .  CHAM .  VOL.  =  "  ;  63 ;  "  CU .  I N .  " 

PRINT  " CHG  NT  =" ; PI ; "LBS. " 

PRINT  "TRAVEL  BIST.  = " ; 62 ; " I H . " 

PRINT  "MUZZLE  VELOCITY  = ,! ;  I  NT  <  V/ 1 2  > ;  "  F*'S 11 

PRINT  "MAX.  CHAM.  PrES.  = 11 ;  I  NT  <  B3  >  ;  "  PS  I  AT  11 ;  £4 ;  "  SEC .  OR"  ;  B5; 11 JN,  " 
PRINT  "MAX.  BASE  PRES.  =  " ; I NT <B8) ; "PS I " 

PRINT  "MUZZLE  PRESSURE  =  " ; I  NT  <  29) ; "PS I  AT  H: SI; “SEC" 

PRINT  "MAX.  ACCEL  = " ; I  NT  < A5 ) ; " GS " 

IF  Z<1  THEN  2130 

PRINT  "FRACTION  BURNT  = " ; 2 ; " AT " ; 25; " IN  FROM  MUZZLE" 

GOT G-  2140 

PRINT  "FRACTION  BURNT  = " ; Z 

PRINT  "BETH" ; B 

PRINT 

PRINT  69  s I  1 , T < 1 , 1 ) , T < 2 , 1 > , T < 1 , 2) , T < 2, 2) 

PRINT 

GOTO  2740 

IF  E9= 1  THEN  2240 

PRINT  " PRO J  NT  BIST  CHAM  VOL  CW  PMA XP 

PRINT  "IV  ACCEL  BETA" 

LET  E9= 1 
PRINT 

PR  I  NT  TAB ( 1 > ; G 1 ; TAB  < 8 > ; 62 ; TAB < 1 5 > ; 63 ; TAB < 24) ; P i ; T  AB < 3?  > : 

PR  I  NT  I  NT < B3 > ; T ; B <40 > ; I  NT  < V/ 1 2 > ; TAB  <48  > ; I NT < A5  > ; TAB <  56 > ; B 

PRINT 

GOTO  2740 

LET  U  =  R < J ,  1 )  • 

LET  X  =  R < J , 2 > 

LET  V  =  R  < J , 3 > 

IF  2 >= 1  THEN  2350 
IF  U >  =  W 1  THEN  2400 

LET  2=1 -<  P8-U *  <  < P7-U  >  A2-P0*  < F9  +  U )  -* 2  >  /Q7 
GOTO  2470 
LET  2=1 

IF  26=1  THEN  2470  j 

LET  25=62- X 

LET  26=1 

GOTO  2470 

FOR  K = 1  TO  12 

IF  F  < 1 ,  K  >  *  W 1  > U  THEN  2440 

NEXT  K 

GOTO  2350 

LET  U3=<  IJ/W 1  -F  Cl,  K-l  >  >  ✓  <  F  <  1 ,  K  )~F<  1 ,  K-  1  >  > 

LET  W4=  <PS-U) /P8 

LET  Z= 1 - W  4  *  <  F  <  2 , K-l >  +  <  F  <  2 , K > -  F  <  2 ,  K  -  1 > >*W3>*W2 
LET  T 1 =0 

IF  X  =  0  THEN  2560 

FOR  K=2  TO  10 

IF  X <  =  T < 1 , K >  THEN  2530 

LET  T 1 =T 1  +  < T < 1 , K ) -T < 1 , K- 1 >  > * < T < 2 , K >  +  T < 2 , K-l > > /2 
NEXT  K 

LET  T2  =  T  <  2 ,  K- 1 )  +  <  T  <  2 , K  > -T  < 2 , K- 1 >  >  *  <  T  < 1 , K > -X  >  ^  <  T  < 1 , K  > -T  < 1 , 1 > > 
LET  T 1 =T 1  +  < T2  +  T < 2 ? K~ 1 > > * < X-T < 1 , K- 1 > > /2 
LET  T 1 =  T 1*64 

LET  T  =  <  J  1  +  Q 1  *  2  -  H  1  * V * V - T 1 >,'< J2+Q2*Z> 

LET  V 0  =  6  3  +  6  4  *  X  -  Q  4  *  < i - 2 > - J 5 - Q 5*2 
IF  V 0 < 0  THEN  2600 
GOTO  2620 

PRINT  G3,G4,Q4, J5,Q5,Z 
GOTO  3030 

LET  B1=T*(Q3*Z+J3)/V0 
LET  E 2  =  B 1  H 2 

LET  E0=B2*H3  30 


t 


I 
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2656 
266  0 
2670 
2680 
2690 
2708 
2710 
2720 
2730 
2740 
2750 
2760 
27?U 
2788 
2790 
2800 
2810 
2820 
2830 
2840 
2858 
2860 
2370 
2880 
2890 
2960 
2910 
2920 
2930 
2940 
2950 
2968 
2970 
2  9  8  0 
2998 

3009 

3010 
3020 
3030 
3990 
4000 
40Q5 
4020 
406O 
4070 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4238 
4240 
4250 
426u 
4  270 

4280 

4281 
4290 

4300 

4301 


IF  XOO  THEN  2670 
IF  B1CG9  THEN  2690 
LEI  A1=H4*<B2-T2> 

GOTO  2700 
LET  A  1=8 

LET  R<J, 1>=2*B*<B1+A> 

LET  R < J , 2)= V 
LET  R<Jf3>=ftl 
RETURN 

IF  H=0  THEN  285© 

IF  D5=0  THEN  2770 
NEXT  N5 

IF  D 4 <  = 0  THEN  2790 
NEXT  N4 

IF  D 3  =  O  THEN  2810 
NEXT  N3 

IF  D2=e  THEN  2830 
NEXT  N2 

IF  D  1  35 0  THEN  2850 
NEXT  Ml 
PAUSE 
DAT A  15 
DATA  185.84,360 
DATA  3,. 00005 
HATH  3400, 3,7.21 
DATA  2500 
DATA  .02 

DATA  1 . 152E+06, 1 . 25,2000,30 
DATA  8.8 

DATA  3445800, 1 . 269,2065,32.81 , . 056- 
DATA  . 1425, . 35, .G132,7 
DATA  "EX 6  5 8 6  " 

DATA  6500, . 000856673, 200,52388,8 
DATA  2,8, 1800,508, 1800 

DATA  1,1,1  ..086, 8. 967,  !  .  024, 8. 874,  1 . 85.5, 0.735 
DATA  1.101,8, 569, 1 . 162,0.395, 1 . 24,8.233 
DATA  1 . 292,8. 191,1. 347,8. 104,1. 489, 0,049 
DATA  1 . 48,0.814, 1 . 557,0 
STOP 

SUE  E j  <  #  1 , #2 , Rho i , V i , T i ,  3H  , T erop i , Ada, Xm . D < C 1 > 

OPTION  EASE  1 

OVERLAP 

Eoa=Eob=0 

X  t  =Xm 

Kount =  T i /Ut 
V i = V i /Xm 
Rh  =  1  /'Rho  i 
T  p  =  8 
G  m  1  = .  2  5 
C  m = . 34 

I  NRIJT  "  Cm ,  Ada,  Grn  1  "  ,  Cm  ,  Ada,  Grn  1 
A 110: K o u  n t  = K o u  n t  + 1 
T  p  =  T  p  +  B  t 
T=T i +  Tp 
D •=■  m  =  1  +  Vi  *Tp 

V  =  V  i  /  li  e  m 
Rho = Rho  i  /Bern 

Temp=T emp i  <Rh~Ada>/< 1 /Rho- Ada) ) A Q m 1 

V  =  V*Xrn 

PRINT  #1 ; V,Xt , Rho, Temp 
IF  Kount  MOD  400  THEN  AllU 
IF  Eoa  THEN  A40 
IF  V*Rho >  = 1 00  THEN  A118 
Eoa=  1 
Dt  =4*D-t 

PR IHT  #1 ; 1E51 ,©,0,0  31 


CO  CO  CO 


4305  GOTO  ft  110 

4310  ft40:IF  Rho>=. 0O050  THEN  ft  1 10 
4315  PRINT  #1 ; 1E51 , 0, 0, O 
4320  X irt 
4330  Dt=Dt--4 

40  PRINT  #2;  Rhoi  ,Vi  ,Ti  ,-Bv,  Tempi  ,  ftda,  Xm ,  I' ,  C 1 
60  SUBEND 
70  END 


t 
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"NOZZLE" 


Input  data 
of  lines  20  and 

Parameter 

F0 

Beta,  PI 

C 

Rhop 

Ada 

TO 

Gamma 

Force 

C0,C2,C2P 

VO 

Ast 

Pstart 

G6 


is  lines  110  and  120  corresponds  to  the  READ  statements 
30.  Input  dat  is  as  follows: 

Data 

Initial  web  fraction  remainging.  (0<F0<1) 

If  F0  is  too  small  the  initial  pressure  will 
exceed  the  intended  nozzle  start  pressure. 

The  assumed  burning  rate,  in  inches  per  second, 
is  Beta  (P+P1)/W  where  P  is  the  chamber  pressure 
in  psi  and  W  is  the  web  size  (inches) 

Charge  weight  (pounds) 

Propellant  density  (lb/cu.in.) 

Propellant  covolume  (cu. in/lb.) 

Adiabatic  flame  temperature  (°K) 

Ratio  of  specific  heats 

Impetus  (in  lbf/lbm) 

These  are  the  kp,  k2,  and  -V-2  of  NDRC  A142  of 
5  Feb  1943;  they  are  used  to  help  calculate  the 
amount  of  propellant  that  has  been  burnt. 

Chamber  volume  (cu.in.) 

Area  of  the  erosion  nozzle  (sq.in.) 

Nozzle  start  pressure  (psi) 

Integration  step  size  (sec) 
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10  OVERLAP 

2  0  R  E  A  D  F  0 ,  B  e  t  a  ,  p  1 ,  C ,  R  h  op,  Ada,  T  0 ,  G  am  rn  a,Fo  r*  ce,  CO,  C  2 ,  C  2  p 

SO  READ  V0 , Ast , Pst  art , G6 
40  S i =S=0 

50  2eta=<2/'<Gamma+l  )  <  1  /'(Gamma*  1  )  > 

Go  Ep-< 1  -  2  *  Z  e  t a > / G am m a 

?0  Phelp=(2/’(G am m a+ 1  > > A < G am m a-- < G am m a- 1 > ) 

80  0  8  -  G  6  -■  2 
90  Break  =0 

1 08  I  MAGE  BB . DDDD , DD . DDL , DDDD . DDDD , DDDDDD , BDDDDDDDD ,  DDL .  DDDDD ,  8D 
j 1 0  IlRT fl  ®  •  500 »  ■  005280 , 6500 ,  .  300 ,  .  057 ,  32 .81, 2065 ,  1 . 269 , 3445000 ,  .  8502 

120  DATA  19. 2, . 098, 41000, . 000050 

130  R  <  1  1  )  =  F  O 

140  R<1,2>=0 

150  Kount=0 

160  CREATE  “DATA", 126 

170  CREATE  "  DATA 1  ",  3 

180  ASSIGN  #1  TO  "DATA" 

190  ASSIGN  #2  TO  "BATR1 “ 

200  BUFFER  #1 

21G  FOR  1=1  TO  2 

220  LET  R  <  2 , I)=R< 1,1) 

230  NEXT  I 

240  LET  J=2 

250  GOSUB  660 

260  LET  S=S+Q8 

270  LET  J=3 

280  FOR  1=1  TO  2 

290  LET  R(2, I )=RC2, I )*G6 

300  LET  R(3, I )=R(1 , I )+R<2, I )/2 

310  NEXT  I 

320  GOSUB  660 

330  LET  J  =  4 

340  FOR  1=1  TO  2 

350  LET  R<3, I )=R<3, I  )*G 6 

360  LET  R<4,  I  )=R<  1 ,  I  )+R<3,  I  >/'2 

370  NEXT  I 

380  GOSUB  660  •  1 

390  LET  J  =  5 

400  LET  S=S+Q8 

410  FOR  1=1  TO  2 

420  LET  R < 4 , I  )  =  R  <  4 , I)*G6 

430  LET  R < 5 , I ) =  R < 1 , I ) +  R ( 4 ,  I ) 

440  NEXT  I 

450  GOSUB  660 

460  FOR  1=1  TO  2 

470  LET  R ( 5 ,  I  )  =  R  <  5  , 1 ) * G 6 

480  LET  R  <1 ,  I  >  =R  <  1  ,  I  )  h-  <  R  <  2 ,  I  )  +R  <  3 ,  I  )  *2  +  2*R  <  4  ,  I )  +R  < 5 ,  I  >  )  /6 
490  LET  R  <  6 ,  I >=RC 1 ,  I) 

500  NEXT  I 

510  PRINT  USING  100; S, 2, Si , T, P, Rho, Ust 

520  IF  U s X  =  0  THEN  210 

5  3  8  K  o  u  n  t,  =  K  o  u  n  t  + 1 

540  PR  I  NT  #  1  ;  List ,  1 ,  Rhost ,  Tst 

545  IF  Kount  MOD  400  THEN  210 

550  IF  NOT  Break  OR  < Rho >.8035)  THEN  210 

570  IF  SI  THEN  SI  a 

580  Sl=l 

590  PRINT  #1; 1E51, 0,8,0 
60O  G6=4"*G6 

605  08=4*08 

610  GOTO  210 

615  SI  a:  IF  Rho >.001  THEN  210 

620  PRINT  #1; 1E51, 0,0,0 

630  PRINT  # 2 J  0 , 0 , 0 ,  G 6 z' 4 ,  O ,  A d a ,  1,  1 . 25*SQR<  Ast  )  ,  1 
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1329, .367 


I 


640  L  0 fl D  " n e w n  s  n “ 

650  STOP 
660  F=R<J,1) 

670  Si=R<J?2> 

68G  2=< l-F>*<C0+C2p*F) 


600  IF  2>1  THEM  2=1 

70O  IF  F >0  THEN  Z=< 1-F>*<C0tC2*F> 

710  IF  F<-.5  THEN  2=1 
720  Zhe  1  p=V0-C*  <  1  -2  >  ^Rhop 
730  Vpres*Zhel  p-<C*Z~Si  !> * fl d a 
749  Vth=2hel p-C*2*flda 

7 bo  T=T0  *  <  Vt  h*  (  2^8  i  /C )  /(, Vpres*Z)  )  <”•(  Gamma- 1  > 
760  P  =  F  o  r  c  e  *  <  C  *  2  --  S  i  >  *T2  <  T0*Vpres  > 

770  Break  = Break  OR  <P>Pstart> 

780  Rho=<C*Z-Si )/2help 


7  9  O  s  i  g  (h  a-  R  h  o  •y  ( 1  —  fl  cj  a«  R  h  o  ) 


SO0  R  <  J ,  1 >  =-Bet  a*  < P+P 1 ) 

810  IF  NOT  Break  THEN  RETURN 

820  He  1 per  =  Ep*flda*Si gma 

830  S  i  gma=.t  =2et-  a*S  i  gma*  C 1  +  He  1  per > 

840  RKost=Si  grnast./<  l+flda*Si  gmast ) 

850  T  s  t  =  2  *  T  *  <  1  +  C  G  am  m  a- 1  >  *  H  e  1  p  e  r  >  /  C  G  am  m  a+ 1 ) 

8t'0  Pet  =Phe  1  p*P*  C 1  +  G  amnia*  He  1  pep ) 

870  U  s  t  =  S  Q  R  ( 3  8  6  *  G  a  m  rn  a  *  P  s  t-  *  (  1  +  R  d  a  *  R  h  o  3 1-  > --S^Rhost  ) 
880  R< J, 2>=fist*Us.t*Rhost- 
890  RETURN 
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"newnsn" 


The  data  listed  in  this  program  should  not  normally  be  changed. 
Input  parameters  are  requested  as  indicated  below.  Default  values 
are  as  indicated. 


Message 


Keyboard  Entry  . 


0  FOR  GRAPHS,  1  FOR  NO  GRAPHS 
TIMM,  TEMPM 

ENTER  LABEL 
GUN  TEMPERATURE? 

ALCP 

XINT,  TINT 


NRDS,  NLA,  IPROPT,  LHELP 


INSIDE  DIAMETER, 
OUTSIDE  DIAMETER 


0  or  1 

Maximum  time  (sec)  and  temperature 
(°K)  to  be  shown  on  graph 
An  alphanumeric  label  for  the  graph 
Temperature  of  gun  steel 
(default  is  300° K) 

(default  is  0.4/(14.8  +  4*logio  D)^, 

D  is  the  gun  caliber  in  centimeters) 
Distance  from  idealized  breech  (inches) 
for  these  calculations,  time  to  stop 
calculations  or  go  to  next  round  (sec) 
(default  =  projectile  base  or  nozzle 
throat,  D/100) 

Number  of  rounds  in  this  burst,  number 
of  mathematical  layers,  print  option 
(0  for  full  print,  1  for  1/4  as  much, 

2  for  very  little  printing) ,  mathematical 
layer  to  use  for  bore  surface  between 
rounds  (default:  1,  7,  0,  4) 

Inside  and  outside  diameter  of  the 
gun  (inches) 
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1  OPTION  BASF  1 

5  DIM  Label TC 40] 

6  I  N  p  U  T  "  0  F  0  R  G  R  ft  P  H  S ,  }  F  0  R  H  0  G  R  ft  P  M  S  “  ,  E  r  i  q 

7  IF  Eng  THEN  35b 

10  PLOTTER  IS  "GRAPHICS" 

11  PEN  2 

12  DEG 

2  O  I  M  PUT  11  T  I M  M ,  T  E  M  P  M  ",Ti  m  m ,  T  $  m  p  n> 

25  L INPUT  “ENTER  LABEL “ , Labe  1  $ 

2  b  Label =  L  E H ( L  ab  e 1  $ ) 

30  LOCATE 

4  0  SCALE  -  T  i  m  m  /  1 0 ,  T  i  m  rn  ,  2  0  O ,  T  e  m  p  m  + 1 0  0 

45  CL  IR.  -T  i  mm/ 1 Q 0 ,  T  i  mm ,  290 ,  T emprn 

5 0  ft  X E  S  ,001, 1 0  0 , 0 , 3 ® 0 

55  UNCLIP 

60  C SIZE  3 

70  LDIR  0 

80  LORG  5 

30  FOR  I = . 0O2  TQ  T \ mm- , OO i  STEP  , 002 

100  MOVE  1,278 

110  LABEL  USING  "HD"; 1**000 

12G  NEXT  I 

121  LDIR  90 

130  LORG  5 

149  F 0 R  1=4 O 8  T 0  T e mp m  S T E P  200 
145  MOVE  - . 030*T i mm , I 
15G  LABEL  USING  “ BPDD 11  \  I 
170  NEXT  I 

1  U  0  L  0  R  G  5 

1 9  0  M  0  V  E  -  f  0  6  *  T  i  m  rn ,  (,  3 0 8  +  T  e  m  p  m )  /  2 

2 0 0  C  SIZE  4 

201  LDIR  90 

218  LABEL  USING  “ 20ft" ; M TEMPERATURE  (DEG.  K > " 

228  LDIR  0 

2  3  8  M  0  V  E  T  i  m  m  /  2  ,  2  3  0 

240  LABEL  USING  " 1 9A" ; "TIME  <M ILL I SECONDS > " 

2  5  O  M  0  V  E  T  i  m  m  /  2 ,  T  e  m  p  m  +  5  8 
260  C SIZE  6 ,  . 5 

270  LDIR  O 
290  LORG  5 
295  PEN  1 

308  LABEL  USING  ,,30A“  ;  Label  $ 

310  PEN  3 
320  FRAME 
338  PEN  4 
348  MOVE  0,300 
355  OVERLAP 

3 61  D I M  Ac  <  9 9  > , ft  1 ( 9 9  > 

37G  D I M  I st  he <  ?  > , Th  < 99 ,4,7), Dr < 99 , 7  > , I st a< 7), I *t o 

,Ak<40),Hold<7>,Kno<7> 

371  DIM  ft h ( 4 ) 

372  MAT  READ  Ah 

373  DATA  . 25, 1 , . 5, . 333333333333 

3  74  I p r  o  p t  =  1 

3  8  O  READ  T  h  t ,  E  o  p  r ,  L  h  e  1  p ,  T  m  a  x ,  ft o  r  d  f 
381  MAT  Hoi d=<Tht ) 

390  DATA  308, 1 , 4 , 6, i 
40G  ASSIGN  #8  TO  “DATA" 

405  BUFFER  #8 

410  ASSIGN  #7  TO  “ DATA 1 “ 

411  Eof=0 

420  READ  #7 ; Rhoi , wi , Ti , Du, Tempi , Ada, Xm, D, Ki nt 

421  Xi=Xint 

430  INPUT  “GUN  TEMPERATURE^" , Tht 
448  MAT  Th“ <Tht ) 

510  MAT  READ  I st a, I st o , Jp , I st he 


J  p  <  7 ) , D  t( 7  > , D  s  <  7  > , ft  1 ph<4@) 
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520  DRTR  1,5, 7,8, 9, 9, 9,7, 11, *13, 15, 15 
530  READ  Ere  ad ,  Er ,  ErO ,  Enew ,  Eos 
540  HATH  1 , 1 , 1 , u, 1 
550  E  r h  =  u 


"560  H2-1E50 
570  Erh=0 


, 15, 15,4 


4,4, 4, 4,4, 4, 1,5, 7,8,9 
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580  R  E  ft  D  S ,  I  r  d  s ,  N  r  d  s ,  K  o  s  hi 

590  DRTR  @,1,1,0 

6O0  fi  1  c  p  =  -  4 <14.8  +  4*LG T ( D >  )  A 2 

610  INPUT  11 ALCP? 11 ,  R 1  c p 

620  INPUT  "HINT,  TINT?1' ,  Hint,  Tint 

630  INPUT  “  NRDS,  NLh  ,  IPROPT,  LHELP?*1 ,  Nrds,  HI  a,  Ipropt,Lhelp 

6 49  T=Q=Tcal=0  f 

73@  L  i  ni  =  Ln*i  i  n=  1 

740  T i nt  ®B*100  * 

750  N  l  a-7 

760  MAT  READ  Ak,A1ph 

770  DRTR  0,0, 0,0,0, . 10$, . 10S, . 107, . 106, . 105, . 103, . 101 , .098, .096, .095, .683, • I 
.073, .073, .068, .062, .053, .@6, .063, .089, . 107, . 122, . 130, . 135, . 140, . 145, . 1 45 
780  DRTR  0,0, 0,0, 0,0, 0,0 

790  DRTR  0, 0,0, O, u, . 12i , . 122, . 1 15, . 109, . 103, .096, .09, .084, .078, .072, .066, .061, 

055, .049,  .046, .045,  .@48,  .@56, .07, .039, . 113, .  123,  .  135, . 148,  .  145,  . 150, .  150 

800  DRTR  0,0, 0,0,0, 0,0,0 

320  Alph«>=0 

338  FOR  k  1  =  1  TO  40 

840  R 1  phwsMAX(fl 1  phfrt ,  R 1  ph<  k  1  >  > 

850  NEXT  K 1 

860  Tdesir*T1nt 

87G  Rlce-422. 5*R1 cp 

880  Hspace=SQR<Du*ftl phm*3 ) 

890  Const  stDu/<Hspac  e^Hspac  e  > 

908  FOP  K 1 = 1  TO  40 

918  A 1 ph<K 1 > ^Ccnst  *A 1 ph  C  K 1 > 

920  Rk O  l )= Const *Rk <kl > 

930  NEXT  K1 

94G  R1  c  p*A  1  c  e*Hspa»:  e*Const 

950  PRINT  "THE  SPRCE  I NTERVAL  IS",Hspac€ 

960  INPUT  "INPUT  INSIDE  DIRM,  OUTSIDE  B I  AM ” , R i , Ro 
970  R i =2. 54*R i 
980  Ro  =  2. 54*Ro 
990  Kou=8-Nl a 


1 0 0 0  S  h  =  H  u  /  H  s p  ac  e 
1016  D  t < K  o u >  = D  u *  1 0 0 0 


1020  Ds  (.  Kou )  =Hspac  e 
1030  FOR  Kol-Kou  TO  6 
1040  Dt  Oso 1  +  1 >=4*Bt <Kol > 

1 050  Ds < Ko 1  + 1 >  =2*Ds  < Ko 1 ) 

1060  NEXT  Ko 1 

1070  I s t  o < 7  >  =  I  NT  <  < Ro-R i )/DsC7>  +  .5> 
1030  REM  --CIRCULAR  CORRECTIONS 
1090  FOR  I  =  K o u  TO  7 
1100  I S  t  Op= I S  t  O  <1  ) 

1UG  FOR  I  i=l  TO  1st  op 
1120  R=R i+Ds<I>*Ii 

1130  DrCIi  ,  I>=Ds<I)/<2*R> 

1131  NEXT  Ii 

1132  NEXT  I 
1140  Eoe=  1 


1145  ftl ce=fl Ice/ 12 
1150  FOR  K 7  = 1  TO  4 
1160  L=7 

1170  Kno< 7>=J  =  K7 
1180  GOSUB  Heater- 
1190  IF  Eoe  THEN  J7 
1200  FOR  K 6  = 1  TO  4 
121G  L=6 


1220  K  n  o  (.  6  >  =  J  =  K  6 
123  U  G O S U B  H  0  a 1 0  r 
1240  IF  Eoe  THEN  J6 
1250  FOR  K5 - 1  TO  4 
1260  L  =  5 

1 2 7  0  K  n  o  <  5  >  -  J  =  K  5 
1280  GOSUB  Heater 
1290  IF  Eoe  THEN  J5 
1300  FOR  K4=l  TO  4 
1310  L~4 

1320  Kno< 4 ) -J -K4 
1330  GOSUB  Heater 

1349  IF  Eoe  THEN  J4 

1350  FOR  K3--1  TO  4 
1360  L=3 

1370  Kno <3 > =J=K3 
1380  GO S U B  H e at e r 
1398  IF  Eoe  THEN  J3 
1 4 0 0  F 0 R  K 2=1  TO  4 
1410  L=2 

1420  Kno < 2 >  = J =K2 
143  0  G  0 o U  £  H  e a t e  r 

1449  IF  Eoe  THEN  J2 

1450  FOR  K 1 = 1  TO  4 
1460  L= 1 

1478  Kno < 1 ) = J=K 1 
148 0  G  O S U B  H  e  a t e  r 
1490  NEXT  K 1 
1508  J2S  NEXT  K2 
1510  J3 : NEXT  K 3 
1520  J4 : NEXT  K4 
1530  J5 : NEXT  K5 
1548  J6: NEXT  K 6 
1  5  5'Q  J  7 :  NEXT  K7 
1560  GOTO  1150 

1570  Heater: IF  Kosh< >0  THEN  2690 

1598  DEF  F N D e  1  <  I  a >  =  ft r  <  I  a , L) ~  <  fik  <  I  s  1  >  *-  fi k  <  I  s  3  >  >  s  <  4 *Rk  <  I  >  >  ; 

1  600  HEF  FNF < H ,  B ,0=01  ( !< > * <R* <  1  -Be )-B~B+C*<  1  tBe >  )+B 
1610  DEF  F N G  <  ft ,  B ,  C ,  E  > = E  *  F  N  F  <  ft ,  B ,  D  >  +  <  1  -  E  >  *  B 
1615  E  o p r =  < I p r  o p t  =  1  )  0 R  < I p r o  p t  =  2  >  R  H  D  < J = 1 > 

1620  Jp<L>=J 
1  6  3  0  R  M  =  0 

1640  Is t art  =  I s t a < L ) 

1650  Is t o p = I s t o < L ) 

1660  Il=Istart+l 

167G  12= I st op- 1 

1680  Jhe=l+< J+2>  MOD  4 

1690  IF  J= 1  THEN  Ho  1 d < L > =Th < I st op , 4 , O 

1700  FOR  1  =  I  start  TO  I  stop 

171G  L  1  =  T  h  <  I  , J h e , L> / 5  Q 

172G  F*  1  =  L  1  MOD  1 

1725  LI =1  NT (LI ) 

1 730  R 1 < I >  =fl 1 ph  <  L 1 ) -P 1  *  <  fll ph ( L 1 > -R 1 ph < L 1 + 1 > > 

1740  flc < I >  =  fi k < L 1 )-Pl *<flk < L 1 >-flk<Ll +1 > ) 

1750  NEXT  I 

1766  IF  L  =  7  THEN  Outside- 
1770  REM  — OUTER  STEP  OF  LRYE*  R 
1780  I s= I st  a(L+ 1 > 

1790  LI =  Th< Is, Jp<L+l > , L+l > / 5 0 
1  8  0  0  K  =  I  s  t  o  p 

1810  De-Dr <K, L)-<fic (I st op-2 > -flk  <X  1 >-Pl >-flk <  L 1 + 1 > >  > / <  4  *  fl c <I*top>> 

1 820  Th  <  I  s  t  op ,  J  ,  L  >  =  F  N  G  <  Th  <  I  st  op-2 , 4 ,  L  >  }  Th<  I  st  op  ,4,0,  Th<  Is,  Jp<  L+  1  >  ,  L+  1  >  ,  fih  <  J  > 

1830  Out  side: FOR  K=I 1  TO  12 

1340  De  =  Dr <K, L) - ( flc <K-1 )-flc <K+i  > <K> > 

1  356  Th  <  K ,  J ,  L  >=FNF  <  Th  <  K- 1 ,  Jhe ,  L  >  ,  Th  <  K ,  Jhe ,  L  >  ,  Th  <  K+l ,  Jhe-  ,L>> 

I860  NEXT  K  39 


1870  IF  L=7  THEN  Th < I st op , J , 7) = <2*Th< I st op- 1 , J , ? > -Th< I st op-2, J , 7 > /2  + 1 . 8 ) / 1 . 566 
1886  REM  --THIS  ASSUMES  K= . 1 1 , H= . 00667, AND  AMBIENT  TEMP  IS  300  K 
1890  IF  L=Lmi n  THEN  A40 
1900  I s= I st  o<L- 1 ) 

1910  U=ThOs-l,  l,L-l>/50 
1920  P 1 =L 1  MOD  1 

1930  De=Dr<  I  start ,  L>-<Ak  (Ll)-Pl  *<flk  <L1  >-flk  <L  1  + 1  >  )-Ac  <  I  start +  1  >  >/<4*Ac  <Istart>> 
1940  K= I start 

1 950  Th  <  I  s  t  art. ,  J ,  L  >  =FNF  ( Ho  1  d<  L- 1)  ,  Th<  1st  art ,  Jhe  ,L),  Th  <  I  st  art  + 1 ,  Jhe  ,  L  >  > 

I960  IF  Eopr  THEN  PRINT  USING  1 955; S, Rv, 0, l, Th< I st art , J, L) , Th< 1st  art  + 1 ,  J ,  L ) ,  Th<  I 
start+2,J,L>,Th<Istart+3,J,L) 

1965  IMAGE  DDD2 . D, 8D2, 6D2. DD , 4D , 7D2 . D, 7DZ. D, 7DZ . D, ?D2. D 

1970  GOTO  Helper  _  .... 

1980  A40: S=S+Dt <L) 

1990  IF  Eread  THEN  A50  .  " 

208G  IF  Eos  AND  <Th<  1 ,  J,  LX  1 . 5*Th<5,  J,  L)  )  THEN  A120 
2010  A45:  ThC  1 ,  J,  L>-<2*Th<2,  J,  L>-Th<3,  J,  L>/2>/l  .  5 
2020  Enew=<S>=Tdesi r) 

2.030  IF  Eopr  THEN  PRINT  USING  1965; S, Rv, Sq, L,  Th< Istart , J , L> , Th< 1st  art  + 1 , J, L) , Th< 
Istart+2, J, L>, Th< Istart+3, J,L> 

2040  GOTO  Helper 

2045  A50:  M=M+ 1 

2046  READ  #8 ; V , X , Rho , Tgas 

2047  X=X+X i 

2050  IF  V  > 1 E  4  9  THEN  A70 
2060  A60: Vgas=0 

2070  IF  X>Xi nt  THEN  Vgas=V*X i  r.t /X 
cd O 8 0  HI =A 1 c e*Rho*Vgas 
2090  Z=H1*H2/<H1+H2> 

2100  Coef=Z*Sh/flc < 1 > 

2110  Th< 1 , J, L)«<2*Th<2, J, L>-Th<3, J, L>/2+Coef*Tgas>/< 1 . 5+Coef > 

2111  R‘,.'  =  2* < Tgas-Th <  1 ,  J,  L  )  > 

2112  IF  Eng  THEN  2120 

2113  IF  S>Timm*1800  THEN  Eng-1 

2119  DRAW  8/1000, Th< Istart, J,L> 

2120  Q =Q+<Rv+R h o v > * D t <L>/2000 

2130  IF  Eopr  THEN  PRINT  USING  1965; S, Rv , Q, L, Th< I st art , J , L > , Th< 1st  art  + 1 , J , L )  ,  Th< I 


£•  t  ar  t  +  2 ,  J ,  L  >  ,  T  h  (1st  ar  t  + 

3, J,L> 

2140 

R  fa  o  v  —  R  v 

2150 

GOTO  H  e* 1 p e r 

2160 

R70: IF  E r  THEN  R80 

2170 

Erea d=0 

2180 

GOTO  fl  6  0 

2190 

fl80 : Er=0 

22O0 

REM - STEP  SIZE 

INCREASE 

2210 

Sh-2*Sh 

2220 

R100:S=S-Dt<L> 

2230 

Lrn  i  n=L  i  ni=Mrri=L+ 1 

2240 

I i =Ista<Mm)-l 

2250 

•J.j  =  K  n  o  <  M  m  ) 

2260 

FOR  1=1  TO  Ii 

2270 

Th  < I , J  j ,  Mrn  >  =  T  h  (  2*  I 

-1,4, Mm- 1 ) 

2275 

NEXT  I 

2280 

I  s  t  a  ( M  rn )  =  1 

2290 

Kosh=3 

2300 

Kount =6 

2310 

GOTO  2366 

2320 

H120:  IF  J<  >1  THEN 

H45 

2330 

Kount = Kount + 1 

2340 

EO'S=  <L<Lhe  1  p) 

2350 

ON  Kount  GOTO  ft  45, 

fl45 , ft  1 08 

2360 

He  1  per : E  o  e  =  <  L  m  i n  > = 

L  > 

2370 

I  st  art  =  1st  a<  L  > 

2380 

I s t o  p=Isto<L) 

2390 

T m a x  =  M fl X ( T rri ax,  T h (  1 

i  J ,  1  >  > 
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2408  IF  Enew=0  THEN  RETURN 

2418  PRINT  ”  Round  no";  Irds;  “ T I  HE  =  M ;  S  s  1 8  0  0 ;  “  Sec  ends  “  ;  U I N  <  1  >  ;  11  Max;  t*R,  pn;  INT<Tma*> 
"Bore  surface  t  ernperat  ure " ;  Th<  1 ,  J ,  Lhel  p>  ;  “Total  heat  i  nput  "  ;  ,  1*INT<  10*Q> 

2428  Q  =  0 
2430  Enew=0 

2440  IF  Nrds.<  =  Irds  THEN  CftLL  Wai  t <NrcJs • Tdes i r , T i nt , Ret urn > 

2450  IF  Return  THEN  RETURN 

24 EG  Tmax=Q 

2478  E o s  =  0 

243  0  I r d  s  =  I r d  s  + 1 

2490  Er ead=Er=ErQ=Tr ue 

2580  T  d  *  s  i  r=Tdesi  r+T  i  nt 

2510  Sh-Sh/2 

2520  Lh=,L-l 

2530  FOR  1=1  TO  Lh 

2  5  4  0  J  p  <  I  >  =  4 

2558  NEXT  I 

2560  I=J 

2578  Mmm=L 

2575  E30 : Lh=M- 1 

2530  I  s  t  a  C  M  rn  m  >  =  1  s  t  h  e  <  M  m  ro  > 

2  598  I  s  =  I  s  t  a  <  M  m  rn  > 

2608  FOR  Ka= 1  TO  Is  * 

2618  ThC2*K-l ,  4,  Lh)=Th<2*K,  1 ,  L  h  >  =  (.  T  h  <  K ,  I ,  M>+Th<K+l ,  l ,  M>  )/2 
262  8  I s t o p  = I s t o < L  h ) 

2630  Ho  1 d <  Lh ) =Th < I s t op , 4 , Lh > 

2640  Mmrn=Lh 
2650  1=4 

2660  IF  Lh>=2  THEN  E38 
2678  Lrn  i  n=L  i  rn=  1 
2630  RETURN 
2  6  9 8  K o s h = K o s  h- 1 
2700  E o e  = 1 
271G  RETURN 

2  7  2  O  S  U  E  W  a  i  t N  r  d  s  ,  T  d  e  S  i  r ,  T  i  n  t ,  R  e  t  u  r  n  > 

2  7  3  8  R  E  ft  Ii  N  a  dd ,  T  a  d  d 
2740  BBT  ft  8,0 

2758  IF  Nadd=G  THEN  BIG  ; 

2760  N rds=Hadd+N rds 
2778  Tint=Tadd 
2788  R  e  t  u  r  n  =  0 
2798  SUBEXIT 

2808  B 1 0 : IF  Tadd=0  THEN  STOP 
28 10  T  des i r  =  Tdes i r+T  add 
2828  Ret urn= 1 
2838  SUBENB 
2848  END 
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